Excess current in superconducting Sr 2 Ru04 
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We present results from point-contact measurements on Sr2Ru04 that show a linear dependence 
of the excess current as a function of temperature and applied magnetic fields over a surprisingly 
wide range of the phase diagram. We propose an explanation of this finding in terms of a p-wave 
triplet-pairing state with coupling to a low-energy fluctuation mode. Within this model we obtain a 
quantitative description of the temperature dependence of the excess current. The impact of surface 
effects on order parameter and excess current is addressed. 

PACS numbers: 74.70.Pq,73.40.Jn,74.20.Rp,74.20.-z 



Introduction. — The discovery of superconductivity be- 
low T c ~ 1.5 K in Sr2RuC>4 has quickly triggered a large 
amount of interest because of the unconventional proper- 
ties pj and the initially proposed analogy to 3 He. The 
enhanced specific heat, magnetic susceptibility, and elec- 
tronic mass indicate the presence of significant correla- 
tions 0,0, IHIll- F° r a more detailed overview see Refs.rn 
and0. The exact symmetry of the superconducting order 
parameter (OP) LZj . I 8tl9UlCl.il ill 121 and notably the pair- 
ing mechanism |13l Il4l Il5l llo . Il7| are still controversial. 
Andreev spectra are sensitive to the order-parameter 



spectra 

symmetry [Tl|,ll8( and are thus an adequate experimental 
probe to yield clarifying information. The shape of the 
spectra previously obtained from point-contact measure- 
ments in superconducting Sr2Ru04 where satisfactorily 
reproduced by an analysis of a p-wave pairing state with 
OP d(k) = z{k x ±iky) d9j. 

It has been shown that the excess current in s-wave su- 
perconductors is proportional to the superconducting gap 
[20| and consequently contains further information on 
the superconducting state. The present paper discusses 
excess-current measurements in Sr2Ru04 in greater de- 
tail. In particular, we present experimental results for the 
excess current in applied magnetic fields and find a sys- 
tematic linear behavior as a function of field over a sur- 
prisingly wide range. This finding is compared with the 
experimental excess-current data from Ref. (squares 
in Fig. which exhibit a strikingly linear tempera- 
ture dependence as well. We show that these two find- 
ings imply a well defined functional relationship between 
the excess current and the OP, and discuss the result- 
ing impli cations in the framework of the p-wave picture 
[71 IllL llil extended to include effects of low-energy fluctu- 
ations |2J, |zj . Our measurements also suggest the pres- 
ence of a normal-conducting surface layer in Sr2Ru04. 
We model such a layer by an enhanced scattering rate 
near the surface and obtain qualitative agreement with 
the experimental point-contact spectra within the p-wave 
picture. 
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FIG. 1: Temperature dependence of the normalized excess 
current across a point contact in Sr^RuCU. Experimental re- 
sults (squares) are taken from Ref.llflt Open symbols show the 
results of our calculation for the excess current from a p-wave 
analysis, without (diamonds) and with (circles) the effects 
of an inelastic scattering channel Fi n (T). The dashed thick 
curve illustrates the scaling relation iexc,2(T) oc AB U ik,2(r)^ 
for the inelastic scattering model. 



Experiment. — Our measurements were performed on 
two single crystals grown both by a floating zone tech- 
nique in different groups. T c was obtained via bulk resis- 
tivity measurements. One crystal, labeled #5 (Ref. EsTl 
shows a midpoint transition temperature T c 50% = 1.02 K 
with a transition width AT C 90% ~ 10% = 0.035 K and 
the other, #C85B5 (Ref. H has T c 50% = 1.54 K and 
^j.90%-10% _ qj 5 k Heterocontacts between super- 
conducting Sr 2 Ru04 and a sharpened Pt needle as a 
counter-electrode were realized inside the mixing cham- 
ber of a 3 He/ 4 He dilution refrigerator. The differential 
resistance dV/dl vs. V was recorded by a standard lock- 
in technique. The differential conductance, dl/dV, is 
obtained by numerical inversion of the measured dV/dl 
data. Measurements were performed in different configu- 
rations with respect to the predominant current injection 
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FIG. 2: Field dependence of the normalized excess current 
across several point contacts in S^RuCU. The magnetic field 
H is aligned almost parallel to the c-axis, and the current 
accross the point contact is applied in the afc-plane. Each 
symbol represents one point contact on one of the two studied 
samples. The full line is a guide to the eye. For explanation of 
-fexc(# = 0) and H Rt (I cxc = 0) see text. The inset shows for 
one point contact typical dl /dV curves from which the excess 
current was determined as a function of magnetic field. 



relative to the crystallographic axis of Sr?Ru04, the ap- 
plied magnetic field, and the surface treatment |25j . Here 
we focus on results obtained for contacts with j\\ab and 
H\\c within an accuracy of about 5-10°. 

Linear field dependence. — We focus on high transmis- 
sion contacts, which exhibit a double- minimum structure 
in the differential resistance, i.e., a double-maximum in 
the differential conductance dl/dV vs. V (see inset of Fig. 
EJ). In this metallic regime, in contrast to the tunneling 
limit, excess current via the mechanism of Andreev reflec- 
tion occurs. Figure shows the observed magnetic field 
dependence of the normalized excess current (symbols) 
across several point contacts in S^RuC^ measured at 
temperatures 0.04 < T/T c < 0.20, together with a linear 
guide to the eye (line). Each symbol represents one point 
contact either on sample #5 (open symbols) or #C85B5 
(filled symbols). The excess current was determined by 
numerical integration of dl/dV vs. V after subtraction 
of a normal-conducting background. The normalization 
values Icxc(H = 0) and H m (I exc = 0) have been deter- 
mined from linear regression of the I exc (H) vs. H data 
for each point contact separately. The inset shows typical 
dl/dV curves from which the excess current was derived. 

The data in Fig. (2J clearly are consistent with a linear 
dependence of the excess current on the applied magnetic 
field. Note that the result is quite universal since it is 
found in both samples in spite of their different value of 
T c . As will be outlined now, the observed equivalence 
of the linearity in the field and temperature dependence 
implies a well defined functional dependence of the excess 



current J exc on the superconducting gap A. 

Scaling relation for 7 exc . — Consider the modulus of 
the superconducting OP near the field dependent critical 
temperature T C (H) as a function of the reduced temper- 
ature t(H) = 1 — T at a given magnetic field: 



A ff (T)| 



t<i 



= A H t v 



(1) 



The mean-field exponent (p-wave approach) is v = 1/2 
while in the recently introduced third-order transition 
picture |26| one has v = 1. The proportionality factor Ah 
depends on the magnetic field, H . In order to find the 
resulting field dependence of the OP modulus consider 
the phenomenological interpolation formula 
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determining the upper critical magnetic field, H c2 , as 
a function of temperature. Eq. J5J with ^lqH c2 (T — 
0) = 1.5 T and T C (H = 0) = 1.5 K reproduces the 
experimental data 27j satisfactorily. The inverse of 
Eq. (J2J determines T C (H). Defining the reduced field 
h(T) = 1 — H at a given temperature and expanding 
Eq. J5J for t <C 1 one finds the relation 
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between the reduced temperature and the reduced field. 
Consequently the reduced field dependence of the gap at 
a given temperature is 



A T (10 



h<l 



A 



H 



[H C 2(T)/2Y 



(4) 



For Ah — constant the pre-factor of the right hand side 
of Eq. Q implies anomalies for low temperatures near 
the critical field, notably lim^^o [limT^o At (H )J_^ 0. 
Since At(^) and J OX c are closely related (Ref. |20j and 
below), the observed linearity in Fig. [21 requires that the 
divergence is compensated by the pre-factor through 
A H ~ [flca(O) - H] v and hence 



constant x A 1 / 1 ' . 



(5) 



Eq. (JSJ marks a central result as it imposes a necessary 
condition on any theoretical approach to the supercon- 
ductivity in Sr2Ru04 in order to satisfy Eqs. and l |4"| 
together with the experimental observations in Figs. ^ 
and |3 In the light of this scaling relation we discuss the 
widely accepted triplet p-wave pairing scenario, which we 
extend to include effects of low- lying bosonic fluctuations. 

P-wave scenario. — The p-wave analysis of the point 
contact spectra applied in Ref. yields the BCS tem- 
perature dependence for the superconducting gap in 
Sr 2 Ru04, shown as AbuIM (thin full line) in Fig.Q] The 
resulting temperature dependence of the excess current, 
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determined from the calculated conductance for Andreev 
type spectra, in the framework of the p-wave analysis is 
also shown in Fig. [T] as I C xc.i (diamonds). The calcula- 
tions are performed for a mean free path of 15 coherence 
lengths (£o = Vf/2irT c ) and for a diffusely scattering sur- 
face modelled as in Ref. It is clear that this model is 
insufficient to describe the experimental data. Neverthe- 
less, it is interesting to note that unlike in the s-wave case 
in unconventional superconductors the excess current is 
not necessarily proportional to the OP; we find near T c 
a temperature variation of the excess current linear in t, 
in contrast to the t 1 ^ 2 variation of the OP. This is be- 
cause impurities and disorder strongly affect the surface 
properties of unconventional superconductors [2g| . 

Pair-breaking by low-frequency bosonic fluctuations. — 
As seen above, the p-wave scenario alone does not ac- 
count for the observed temperature dependence of the 
experimentally obtained I oxc (T). This is true also for 
the overall magnitude of the bulk gap, A(0) = 1.1 meV = 
6 x 1.76/cbT c , extracted from tunneling spectra 0]. To 
reconcile the measured A(0) and I exc (T) with a p-wave 
OP we consider an additional pair-breaking channel. It 
was shown by Millis et. al. [21| that a low-frequency 
bosonic mode at a characteristic frequency oj p described 
by an Einstein spectrum A p (lo) = ^J p uj p S(uj — uj p ) leads 
to a temperature dependent pair-breaking parameter 

T in (T) = ^lj p c p coth(^), (6) 

where g is the coupling-constant appearing in the gap 
equation. The assumptions are that tu p < T c <C cl>e, 
we being the frequency of the pairing mode, and that 
A p (lo) is unaffected by the transition into the super- 
conducting state. We performed calculations using the 
quasiclassical Green's functions technique and included 
the pair-breaking parameter as a self-energy within a 
self-consistent Born approximation, i.e. Ej n (i?, e, T) = 
r in (T) (g(pj, R, £))p f ! where e is the energy of the quasi- 
particles, R the position with respect to the interface, 
and the Fermi momentum; the pj-average is over 
the Fermi surface. The Green's function g(p^,R,e) is 
a functional of the self-energy Y,- m (R, e,T) in the usual 
way. The order parameter profile A(i?, T) near the in- 
terface was then obtained by iterating the weak-coupling 
gap equation and E; n (ii, e,T) until convergence. 

For the excess current this model gives an excellent 
agreement with experimental data, as shown by / C xc.2 
(circles) in Fig. □ for u p = 0.5T C and ^^-J P = 2tt x 
0.25. The almost linear temperature dependence over the 
whole temperature range is reproduced within our model, 
and furthermore, as shown as the dashed thick line in 
Fig. ^ the above introduced scaling relation between the 
calculated 2exc,2(T) and the theoretically obtained order 
parameter AB U ik,2(T') is fulfilled to remarkable accuracy 
with the scaling exponent v = 1/3. 



Another effect of Tin (T) is that the enhanced scattering 
at higher temperatures reduces the observed T c substan- 
tially from its ideal value while the gap at T — > is much 
less affected, giving A(0)/fcBT c -ratios much larger than 
the BCS-ratio 1.76. Our calculations give the correct 
absolute magnitude, A(0) = 5.6Abcs(0). Notably, also 
the functional form of A(T)/A(0) is modified compared 
to the pure p-wave case (see AbuIM, thick line in Fig. 
P). The conductances calculated with the present model 
have the same qualitative features, both for the Andreev 
and the tunnel limit, as those displayed in Ref. 0], and 
can still account for the measured data. 

Normal-state surface layer. — In order to obtain more 
detailed insight about the nature of the pairing state in 
Sr 2 Ru04 it would be instructive to quantify empirically 
the field dependence of An in Eqs. and Q. Un- 
fortunately, obtaining data from the necessary temper- 
ature scans at different fields for a given point contact 
is difficult because of the sensitivity of the large back- 
ground resistivity 0] in the dV/dl data to very small 
changes in the configuration. A possible reason for the 
presence of a large background can be found in a normal- 
state surface layer due to surface reconstructions 29] that 
leads to an additive resistivity in the point contact as 
^measured = #N + °- N i s - Here ct n 1 s is the conductivity of 
the normal-superconductor interface, is the normal 
layer resistivity, and typically cr^i s /i?N ~ 10%. Note 
that the thickness of the normal-state surface layer ap- 
pears to be independent of the sample quality since the 
observed values of 0.5 fl < i?N < 25 f2 vary from point 
contact to point contact but are in the same range for 
both samples |2fij . 

Such a normal-state surface layer has a natural expla- 
nation in a p-wave triplet scenario because the p-wave OP 
is very sensitive to scattering. We assume a region near 
the interface in which scattering is enhanced. In Fig. |3 
we show the self-consistent OP, A(R,T = 0.05T C ), for 
a mean free path of 0.3 coherence lengths in the shaded 
region, and of 15 coherence lengths elsewhere. The bulk 
OP is of the form k x + ik y , and near the surface a sec- 
ondary OP component, k x — ik y , is induced. As can be 
seen from Fig.01 both components are suppressed in the 
surface layer where scattering is enhanced, leading effec- 
tively to a normal-conducting layer near the interface. 

The presence of a normal-state layer affects strongly 
point-contact and tunneling spectra. However, as we 
show in the insets of Fig. 01 the form of the spectra in 
the presence of a normal-state layer is in agreement with 
experiment (c.f. inset in Fig. [21 and Ref. Il9|) . The excess 
current is reduced by such a surface layer (see lower inset 
in Fig-OJ) • Also, the tunneling spectra show a pronounced 
zero-energy anomaly in contrast to the clean surface. The 
temperature dependence of the excess current near T c 
is only weakly affected by a normal-conducting surface 
layer, leaving the results discussed above unaltered. 

Conclusions. — We presented point-contact measure- 
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FIG. 3: Creation of a normal-conducting surface layer in a p- 
wave superconductor due to an increased scattering rate near 
the surface. For comparison, as dashed lines are shown the 
curves assuming a clean surface. The two OP components 
are the bulk k x + ik y order parameter, and the subdominant 
k x — iky order parameter which is stabilized only within a few 
coherence lengths (£o = w //27rT c ) near the surface. Both OP 
components are suppressed in a layer with increased scatter- 
ing, leading effectively to a normal-conducting surface layer. 
The calculations are for T = 0.05T C . The insets show the 
corresponding point-contact spectra (bottom) and tunneling 
spectra (top). 

merits on the unconventional superconductor Sr2Ru04 
as a function of temperature and applied magnetic field. 
The excess current exhibits linear behavior both in tem- 
perature and magnetic field over a large range. Using 
these findings we derive a scaling relation between the 
excess current and the order parameter [Eq. JSJ)]. 

We discuss this result within the theory of p-wave spin- 
triplet superconductivity. We find that that the excess 
current in unconventional superconductors is not neces- 
sarily proportional to the order parameter. In order to 
account for the wide range over which the linear behav- 
ior of the excess current holds experimentally we extend 
the pure p-wave theory to take into account scattering 
between quasiparticles and low-energy bosonic fluctua- 
tions, probably originating from spin fluctuations |2l| . 
The extended theory yields a very good agreement with 
the measured excess current and yields a scaling expo- 
nent v — 1/3. Furthermore it can account for the large 
A(0)/fcgT c -ratios obtained from point-contact measur- 
ments 19] . Finally, we show that surface effects should 
be considered for a satisfactory reproduction of the point- 
contact spectra. 

In closing, we mention that a recent Ginzburg-Landau 
analysis, assuming a third order phase transition induced 
by gapless excitations in the superconducting phase, 
yields the correct temperature dep endence to account for 
the data, at least close to T c |22j, |26( . As shown in Ref. 
8] the p-wave channel of superconductivity may be only 
marginally dominant assuming that pairing in Sr2Ru04 



is mediated by incommensurate spin fluctuations. In this 
case the presence of fluctuations in the OP is not unlikely. 
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